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SUMMARY 


This portion of the report explains the use of the computer program that 
is described in Part 1 of the report. Certain detailed calculation procedures 
are documented. However 9 the formulas of Part 1 are not duplicated here. The 
main purpose is to explain the input to the program so that tha reader is pre- 
pared to use it for his application. 
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I. INTRODUCTION 


An exact general method of solving the Neumann (or Second Boundary 
Value) problem has been applied to the calculation cf low speed flows about 
or within bodies of axially symmetric shape. Solid body, irleu, and purely 
internal flow problems can be solved. The method is capable cf dealing with 
several bodies at once in the presence of one another, and consequently, 
interference problems can be treated with ease. Boundaries need not be 
solid, that is, flows involving area suction can be calculated. Velocities 
can be computed not only for points on the body but for the entire flowfield. 

A surface source distribution is used as the basis for solution. This 
leads to a Fredholm Integral Equation of the second kind, which is solved 
as a set of linear algebraic equations. The source distribution can'be^ taken 
as constant, linearly varying, or quadratically varying across each of the 
elements which form the body. The elements themselves can be flat or Curved 
(parabolas). 

At the present time, it is possible to obtain the solution for the flow 
about any body that has the previously mentioned characteristics and whose 
profile can be defined by no more than 500 coordinate points. 

This method is a generalization of a previous program [1-6]. In 
particular the analysis of part 1 of this report is incorporated into the 
program of [5]. 



II • EQUATIONS AND METHOD OF SOLUTION 


A, Formation of the Elements 


Each body Is described to the present program by the input of an ordered 
sequence of points distributed along that body f s surface* By convention, the 
flowfield of interest is assumed to be "overhead" to an observer who Is 
"walking 11 along the body surface in the order of the input point sequence* 

This convention is necessary in order to establish the direction of the "out- 
ward’ 1 normal vector* 

Surface elements are constructed by connecting consecutive points in the 
sequence by either straight lines or parabolas. (Recall that since the body 
is assumed to be axisyrmnetric , and if we take the x~axis to be the axis of 
symmetry, these lines or arcs should be thought of as "conoids" formed by 
their rotation about the x-axis). In the case of straight (hereafter "flat") 
elements, each element is simply determined by straight lines constructed 
between consecutive points* However, for the case of parabolic arc (hereafter 
"curved") elements, the curvatures may either be Input by the user or deter- 
mined automatically by the program in the following way: the curvature of 

each element (except for the first and the last elements of each body) is the 
geometric mean of the curvatures of the two circles defined by the two end- 
points of each interval and the next sequential point on either side of that 
interval. For the first and last elements, the curvatures are simply those 
of the circles defined by the first three and last three points of the body, 
respectively. The curvature, of each of these circles, is calculated from the 
formulas of Appendix A of Part 1 of this report. 

K = d/r , 


where. 


d = 4[(x 2 -x i ) (y 2 -y 3 ) - ( x 2 -x 3 Hy 2 -y ^ 1 


and 


r = Y(x 2 d-d h ) 2 + (y 2 d ~ d k ) 4 


where. 


^ = 2[r i 2 (y 3 -y 2 ) + ^(y^) + ^(y^)] 


and 


\ = 2[r^(x 2 -x 3 ) + 
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where, 


r j 2 = + y J 2 

r/ = x 2 2 + y 2 2 

r 3 2 = X 3 2 + ^ 


Once the correct geometric means of these curvatures are calculated to 
determine the local element curvatures, the control point locations of each 
element are obtained by approximating the circular arc with a parabola having 
that curvature at its vertex (this is correct through third order in As) • 


B* Calculation of the Influence Coefficient Matrices 

The velocity induced by an element at a control point is computed by 
integrating the basic ring source (or vortex) formulas over the element as 
described in Fart 1. For j i 1 the integrals of equation (9) of Part 1 
are evaluated by Simpson ? s rule in the following manner. 

The integrand is evaluated at an odd number NI of points equally 
spaced in s over the element, where 


NI = Integer Part of j 2 ^8 + 0.9]+ 1 j (1) 

where As is the total arc length of the element and d min is the distance 
between the point where velocity is being evaluated and the centroid or one 
of the endpoints of the elements, whichever is closest. (If NI turns out 
to be less than 3, it is set equal to 3; and if NI turns out tjj be greater 

than 129, it Is set equal to 129*) -Once these NI values of V have been 

determined, the Simpson rule integration formula is applied: 


where 


f V ds = ~ [v x + 4V 9 + 2V 3 + 4V 4 + 2V 5 + ■ 

As^ 

] 


+ 2V + 4V + v 
NI-2 NI-1 NI 


( 2 ) 


V x = vjx^ y i# b(-Sj'), a(-Sj'),J 
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v 2 = V^, y ± , b(-Sj *+6s) , a(-Sj'+6s)J 


etc. 


and where 6s = Sj/(NI-1). 

The second integral in equation (9) of Part 1 is now easily determined 
by multiplying each of the (V , m « 1, NI) values by the local value of s, 
and then apply the Simpson rule formula to the results 



Vs ds = ■y[(V 1 )(-s j ') + 4 (V 2 )C-s 1 , +«8> + 



This procedure is repeated to evaluate the third integral in equation (9) . 


It should be pointed out that equation (9) of Part 1 encounteres diffi- 
culty when the parameter: 


■yStXj-b j) z + aT 1 




(4) 


where e = 0.01, When this happens (termed "small y"), expansions of these 
equations (in this parameter) are used. Their expansions are in reference 3. 

For j = i, use is made of the singular— subelement expansions of 
Appendix B of Part 1. 


C. Matrix Assembly 

The assumption is made that the surface source strength, a, varies 
parabolically between three successive values of a. The derivatives, 
and may then be expressed in terms of the local source strength values 

at the control points of the elements on both sides and then may be replaced 
in equation (10) of Part 1 as follows: 


i 


V * = 
ij ij J 


+5 i®[ A u/„-l <0>+ A 12j°M® + A 13j°M+l f0) ] 
4? if[ A 21j’H-‘ 0,+ A «J ^ >+ A 23 ”M +l W ] 


( 5 ) 




where. 


j + 1 
j - 1 


for the first element of each "body" 
for the last element of each body 
for the other element of each body 

for the first element of each body 
for the last element of each body 
for the other elements of each body 


A 13j 


( g 2 

V %-2 

f Ga 


for the first element of each body 
for the last element of each body 
for the other elements of each body 

for the first element of each body 
for the last element of each body 
for the other elements of each body 

for the first element of each body 
for the last element of each body 
for the other elements of each body 

for the first element of each body 
for the last element of each body 
for the other elements of each body 


for the first element of each body 
for the last element of each body 
for the other elements of each body 



6 




.{ 




L j 

The above A., .*•, L, are functions only of the geometry , and are 
calculated prior io initiating the actual Vjj matrix formulation. The pro- 
cedure for handling the Vij matrix to produce the normal- and tangential-flow 
matrices is identical to that eniployed in the original program (see reference 5). 
The vector velocity is dotted into the local element normal (u^) and tangen- 
tial (&£> directions 



A ij = V ij * Vj “ V ij . * ( ' sin “j C x + cos «j V 

T ij = ' V = \j * (coS “j °x + “j Q y > 


The normal velocity matrix is then used to calculate the source strength values 
for each onset flow as: 




( 9 ) 


7 



~T 


where V» . is the onset normal flow at the i-th control point. For the axisym- 
metric onset flow: 


V, 


°°i 


sin a. 


( 10 ) 


and Vy is determined as described in Section B.l, For the crossflow case as 
explained in reference 2, the normal onset flow truly varies as the cosine of 
the circumferential angle; but the source density, which is truly cosinusoidal 
may be obtained by applying the zero normal flow requirement at 0 = 0. For 
this case, the right-hand side of equation (9) is found from: 


V ooj[ = cos o^. (11) 

Note that the Ay matrix for the crossflow case is different from that of the 
axisymmetric case. 


D. Solution of the Matrix Equations 

The solution of equation (9) for the axisymmetric onset flow and for the 
crossflow cases is obtained by inversion of the A^. matrix, as explained in 
detail in reference 6. This procedure allows multiple "right-hand sides" to 
be computed with very little increase in progr^jn run time (beyond that of a 
single right-hand side) . 


It should be realized that the constant surface vorticity option auto- 
matically generates a number. jof non-uniform axisymmetric onset flows equal to 
the number of "bodies" that have^heen-Jjiput , 


E. Calculation of the Surface Velocities 


Once the source strengths have been determined, the net surface veloci- 
ties are calculated as follows: 


V i 


3=1 


V, .a. + V^. 
Ij j ~i 


( 12 ) 


For axisymmetric flows, the velocity is presented in tangential component 
form. For the crossflow case, however, as discussed in reference 2, the 
velocity is presented in terms of the "fundamental velocity components, T 
and T^,* 1 where 2 



( 13 ) 
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Thus, the velocity tangent to a meridian curve (at any 0 value) is cos 6* 
The "crossflow" velocity (i.e«, the velocity component tangent to the circular 
cross-section) is -T^ sin 0. 






III. INPUT INSTRUCTIONS 


A. Card Column Parameter Locations 


Card 1 - 

- Header Card* 

ccl-60 

Header 

Any alphanumeric run description. 

cc63-68 

Case 

Case number. 

Card 2 - 

Control Flag Card 

ccl 

NB 

Number of bodies (1<NB<9). 

cc2 

NNU 

Number of non-uniform flows (0<NNU<5). 

cc3 

IAXI 

Axisymmetric flow flag. 

cc4 

ICR0SS 

Crossflow flag. 

cc5 

I0FF 

Off-body point input, flag. 

cc6 

I0NLY 

Basic-da ta-only flag. 

cc7 

IELPSE 

Ellipse generator flag (see also Card 5). 

cc8-10 

(blank) 


cell 

IPRTRB 

Perturbation velocities only. 

ccl2 

IP0TNL* 

Solve potential matrix. 

ccl4 

IPTANV 

Prescribed tangential velocity (for the last 
IPTANV bodies) . 

cc!5 

IV0RT 

Strip— ring vorticity flag. 

ccl6 

I0MITA 

Omit axisymmetric uniform flo** T solution. 

ccl7 

I0HITC 

Omit crossflow uniform flox* solution. 

ccl8 

ISURFV 

Surface vorticity (instead of sources) for the final. 
ISURFV bodies. 

ccl9 

IPRSCV 

Prescribed values of the surface vortex strengths 
for the final ISURFV bodies will be input. 

cc20 

IALLV 

All bodies are surface vorticity bodies. 

cc21 

IEXCRS* 

Extra crossflow. 

cc22 

IGENBC 

Generated boundary conditions. 

cc23 

IRNGW 

Ring wing option. 

cc23 

IPNCH 

Punched output. 

cc29, 30 

IUNIT 

Unit number for input coordinates (default = 05) . 

cc31 

IVIJ 

Matrix print flag. 

cc32 

IC0EF 

Matrix-assembly coefficient print flag. 

cc33 

IPRINT 

Very detailed matrix construction print flag. 

cc34 

IRAKF 

Automatic rake generation flag (see also Cards 8 & 9) 

Card 3 - 

Chord/Mach 

number card. 

ccl- 10 

CH0RD 

Reference chord length (default ~ 1.0). 

ccll-20 

XMACH 

Mach number for Goethert correction (0.0 implies 
incompressible) . 

cc21-80 

(blank) 



* Available if and only if N0NEWF * 1, ISIGF * 1 and IGE0MF - 1. 
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Card 4 - Body Control Card 1 of 2 


ccl 

IGE0MF 

0 « curved elements; 1 « flat elements. 

cc2 

ISIGF 

0 « parabolic o; 1 « linear a; 2 « constant a 
(on each element). 

cc3 

ICURVN 

0 « internally calculated element curvatures; 
1 = input curvature (see card 7). 

cc4 

N0NEWF 

0 * use the newest formulae; 1 - use the old formulae 
(implies flat elements and constant a) . 

cc5 

IF0RMT 

Input format flag (see Card Set 6) 

cc6-10 

NN 

Number of defining endpoints for this body. 

ccll-20 

XMULT 

x-multiplier value (default « 1.0). 

cc21~ 30 

YHULT 

y-raultiplier value (default « 1.0), 

cc31-40 

THETA 

Coordinate rotation value (degrees, measured about 
- Z-axis). 

cc41-50 

ADDX 

x-increment (to be applied to all the input coordinate 
for this body). 

cc51-60 

ADDY 

y-increment (to be applied to all the input coordinates 
for this body) . 

Card 5 - 

Body Control Card 2 of 2 

ccl-10 

IBDN 

"Body” number (sequential for bodies, zero for 
off-body points) . 

ccll-20 

ISUBKS 

Subcase flag. 

cc21-30 

NLF 

Non-lifting flag (for combination cases, only) 

cc31-40 

A 

Semi-major axis for ellipse cases) 

cc41-50 

B 

Semi-minor axis for ellipse cases j ^ IELPSE^O. 

Card (Set) 6 - Body 

Definition Cards 

IF 

IF0RMT=O : 

X-coordinates (6F10.5), then 
Y-coordinates (6F10.5) 

IF 

IF0RMT=1: 

X, Y coordinates (2F10.5) (i.e„, one "point-set" 
per card) 

IF 

IF0RMT=2 : 

X, Y coordinates (FlO* 5, 10X,F10. 5) (i.e. , one "point- 
set" per card). 

Card (Set) 7 - Input curvature values (needed only if ICURVN ^ 0), 

(6F10.5) 

CURV(I) ,1= 

=1,NN-1 The curvature values for the NN-1 elements 


which constitute this body. 


Repeat Cards 4-7 a total of (NB+I0FF) times 


Card 8 - Rake Number Card (needed only if IRAKF ^ 0) 

The number of "automatically” generated mass-flow 
rakes. 


ccl-10 


NRAKES 



Card 9 - Rake Definition Card (needed only if IRAKF 4 0) 


ccl-10 

Xll 

cc 11-20 

Ylt 

cc21-30 

X2 1 

cc31-40 

Y2t 

cc41-45 

N 


Coordinates of "start" of the rake. 

Coordinates of the "end" of the rake. 

Number of intervals to be used in the rake (note, 
4SN<200, and N must be an even integer). 


Repeat Card 9 a total of NRAKES times. 


Card 10 - Non-Uniform Flow Control Card (needed only if NNU 4 0) 


ccl~10 NUN 
ccll-20 MSF 

cc21-30 TYPE 
cc31-40 FG 


Flow identification number. 

0 = axisymmetric onset flow; 1 crossflow onset flow, 
2 * both 0 and 1. 

+1.0 = velocity will be input in x,y component 
form; 0 - velocities will be input in normal, 
tangential form; -1 = automatic generation of 
flow due to rotation about the Z-axis (for 
crossflow, only). 

Flow generator constant. 


Card (Set) 11 - Non-Uniform Flow Velocities (needed only if NNU ^ 0) 

(6F10.5) VX(I) or VN(X), I - 1 , total number of control points. 

(6F10.5) VY(I) or VT(I), 1=1, total number of control points. 

Repeat Cards 10 and 11 a total of NNU times. 



Card (Set) 12 — Specified Tangential Flow Velocities (needed only if 
IPTANV 1 0 

(6F10.5) TG(I), 1=1, total number of control points on the last 
IPTANV bodies 

B . Discussion of Parameters 

In all c«*£s, the integer values (indicated by the standard F0RTRAN IV 
naming convention) must be right adjusted within the specified input fields. 

For floating-point type input, the input decimal point will override, the quoted 
F0RMAT specifications. 

Card 1. Self-explanatory. 
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Card 2. Although most of these flags are unchanged from the original program, 
a brief description of each one will be given here for completeness. 

NS The total number of separate "bodies" (1 < NB £ 8) . 

NNU The total number of (user specified} input non-uniform flows 
(0 < NNU £ 5) . 

IAXI Axisymmotrit flow flag (G: no axisymmetric flow, 1: axisymmetric flow) 

ICROSS Crossflow flag (0: no crossflow, 1: crossflow) 

IOFF Off-body point input flag (0: no off-body points will be input, 1: off 
body points will be input) 

IONLY Basic-data-only flag (0: full execution, 1: basic data only) 

IELPSE Ellipse generator flag (see also card 5)(0: standard body defining 

coordinates will be input, 1: coordinates for the ellipse whose proper- 
ties are specified on card 5 are to be created automatically by the 
program) . 

IPRTRB Perturbation velocities only (0: standard total net velocities includ- 
ing the onset flow will be printed, 1: perturbation velocities only will be 

printed) . 

IP0TNL Solve the potential matrix (0: standard solution for the surface 

velocities, 1: solution for the surface velocity potential). Since 
the velocity potential problem for the higher order elements has not 
been coded in the present version, the solution for the potential is 
available only for the case of flat elements and, in particular, 
through the usage of the "old" velocity formulae. These "old" formulae 
have been retained as a subsection of this new, higher order program 
and may be "reached" by setting the following quantities (on card 4) * 
IGE0MF = 1, IS IGF = 2, and N0NEWF = 1; in such a case, the solution 
for the potential is still available. 

IPTANV Prescribed tangential velocity flag (0: standard zero normal velocity 

solution is to be obtained, 1: a user defined set of prescribed tangen- 
tial velocities will be input on the last IPTANV bodies; obviously, 
IPTANV ^ NB.) See the Appendix for further details on the use of 
this flag. 

IV0RT Strip ring vorticity onset flows for each of the bodies (0: no strip 
ring vorticity, 1: automatic generation of the strip ring vorticity 
onset flows) . 

I0MITA Omit the uniform axisymmetric flow solution (0: calculate the standard 
uniform axisymmetric flow solution if IAXI ^ 0, 1: omit this uniform 
axis yuniie trie flow solution even if IAXI £ 0). 

I0MITC Omit the uniform crossflow solution (0: calculate the uniform crossflow 
solution If ICR0SS ^ 0, 1: omit the uniform crossflow solution even if 
IGR0SS ^0) 

ISURFV Surface vorticity for axisymmetric flow (in place of surface sources) 
flag (0: use the standard source distribution on all elements, ¥ 0: 
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IPRSCV 

IALLV 

IEXCRS 

IGENBC 

IRNGW 

IPNCH 

Card 3. 
Card 4. 


use vorticity instead of a source distribution as the singularity on 
the elements of the last ISURFV bodies). Note that if ISURFV ^ 0, then 
both IV0RT and IPTANV must be ^ 0, and then in particular, IPTANV must 
« ISURFV). If we assume that IPTANV « ISURFV = k, and that M is the 
total number of elements on the last k bodies, then the surface vor- 
ticity option causes the source induced .velocity formulae to be used 
on the first N-M elements of each matrix row, and the vortex formulae 
to be used on the remaining M elements of each matrix row. The solu- 
tion for the unknown strengths then proceeds as it would have if only 
IPTANV was nonzero. ; 

Prescribed vorticity flag for axisymmetric flow (0: prescribed vorticity 
values will not be input, 1: prescribed vorticity values (Wj) for the 
last IPRSCV bodies will be input; note that if IPRSCV is ^ 0, then both 
IV0RT and IPTANV must ^ 0, and then, in particular, IPTANV must « IPRSCV). 
If IPRSCV is 0, thea the vortex strengths of the affected elements are 
taken to be 

Total vorticity flag for axisynmetric flow (0: all elements are assumed 
to be source type elements [unless ISURFV 4 0] , 1; all elements are 
assumed to be vortex type elements). Note that if IALLV is £ 0, the 
use of IPTANV is optional, but not mandatory. 

Extra crossflow flag (0: no "extra" crossflow, 1: generate an "extra" 
crossflow [i.e., having a potential which varies as the cosine of twice 
the circumferential angle]). Note that if IEXCRS ^ 0, then N0NEWF - 1, 
ISIGF and IGEOMF must = 1 , since the higher order formulae for this 
kind of velocity potential have not been included in this program. 

Generated boundary condition flag (0: do not generate any additional 
boundary conditions for the crossflow case; 1: generate the onset cross- 
flow due to rotation about an axis normal to the axis of symmetry [see 
also the notes for Card 10]). 

Ring wing option (0:no ring wing option, 1: use the ring wing option 
[see MLC Report J074i/01, April 1970 for further details]). 

Punched card output flag (0: no punched card output, 1: punched card 
output) . 

This card is usually left blank, which results in a default chord 
length of unity, and no Mach number corrections (i.e., incompressible 
results). If a non-zero value of the Mach number is input, the 
program uses the Goethert correction to account for compressibility. 

Usage of the IGEOMF and ISIGF flags permit the user to "turn-off" any 
or all of the higher order element curvature and/or varying source 
density terms. The default values are curved elements with paraboli- 
cally varying source density. If IGE0MF =0, the program will auto- 
matically calculate the local element curvature values by the procedure 
described in Section II .A, unless ICURVN is non-zero, in which case 
the user must supply these curvature values (see Card 7). The value 
for the N0NEWF flag is ordinarily left blank (or zero), even if a flat 
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Card 5 


Card 6 

Card 7 
Card 8 
Card 9 


element, constant source density solution is required* However, since 
certain of the original Doug las -Neumann Axisymraetric Potential Flow 
Program capabilities have not been made available in the higher order 
program (e.g. , calculation of the potential [hence the added mass, etc*) 
the original formulae have been preserved within this version (see 
description of the IP0TNL and IEXCRS flags) . The original capabili- 
ties can be obtained by setting IGE0KF - 1, I3IGF =2, and N0NEWF = 1 
for each of the input bodies, in which case the input instructions of 
reference 5 apply. 

The other parameters on Card 4 are self-explanatory , with the understanding 
that the order of coordinate transformations areas listed on the input. 

The "bodies” are normally loaded prior to any off-body points 
(although the latter require the usage of body control Cards 1 and 
2, also). For this reason, the value of ISDN should be sequentially 
increasing beginning with unity. A non-zero value for IStTBKS means 
that the body definition cards for this body (Card Set 6), (say this 
is the ith body input under this header card) will not be included, 
but that the program is to use the ith set of points that were input 
under the previous header card. Obviously, this capability is useful 
only with "stacked" input cases . 

The three input formats that are available are as shown. The value 
of IF0ROT determines which format is used. The default format is 
the "old" format: x-ccordinates , followed by the y-coordinates . 

These values are to be input only if ICURVN ^ 0. 

This card is needed only if IRAKF J 0. Note that 1 < NRAKES < 20. 

These are tfee rake definition cards, which define the "start" and 
"end" of each mass flow rake. The sign convention that is employed 
is as follows: a positive mass flux means that the flow is from 

left to right to an observer traversing the rake from point (XI, Yl) 
to point (X2, Y2) . For example, for an inlet without a bullet, point 
1 is typically at y = 0.0, and point 2 would have a y-coordinate 
located on the inlet wall. Note that the x-coordinates need not be 
the same, i.e., a tilted rake may be used, if desired. Note also 
that one rake may "overlap" another if so desired, since they are 
each treated independently. It should be pointed out that a Simpson 
Rule integration is performed over the "N" intervals (and therefore 
N must be an even integer) , with the program automatically generating 
N-l "intermediate" points as "off-body points" for each input rake 
(program limit is 500 total off-body points). The velocity values 
at the first and last rake points are obtained by linear extrapolation 
of the two nearest values. In this way, the difficulty associated 
with calculating induced velocities at off -body points which lie very 
near to, or on, the surface is avoided. Typically, values of N between 
10 and 20 appear to be satisfactory for most cases. 
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Card 10 


• Self-explanatory, except that FG, if entered is used in the following 
way to generate the rotation onset flows 

= yi 

V y± = FG-x ± 

v n “ v > : i sin °i“ v yi cos “i 
V T = Vx i COS a i + V yi sin a i* 

Card 11. Self-explanatory. 
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IV. OUTPUT FORMATS 


The main output of the program, which consists of calculated surface 
velocities and/or pressures has not been changed. Its format is that of 
[5]. The surface coordinates with which the calculated velocities are 
associated are the transformed coordinates, i.e., the original input coordi- 
nates altered by uiy specified translations, multiplications, or rotations. 
Similarly, the format of the calculated velocities at off-body points is 
similar to that of [5]. The only change is that output at rake points 
(q.v.) is included. 

The initial output, which details the surface geometry, has been 
changed considerably. The header identifies what kind of solution has been 
computed: (1) curved or flat elements; (2) constant, linear, or parabolic 

source density element curvatures. The header also specifies whether "new" 
or "old" velocity formulae have been used, the latter of which apply only 
to the flat-element constant-source case. In the body of the output, the 
first two columns are the untransformed (input) coordinates, and the third 
and fourth columns are the transformed coordinates, which are the endpoints 
of the surface elements. The fifth and sixth columns are the coordinates 
of the control points of the elements. Element arc lengths are given in 
column 7 and a running total of arc length is displayed in column 8. 

Column 9 lists the differences between average slopes of successive elements, 
i.e., the differences in the slopes of the two straight lines, through the 
respective element endpoints. The actual slope discontinuity at an endpoint 
between two parabolic elements is normally much smaller than the difference 
between the average slopes of the two elements. Moreover, this discontinuity 
approaches zero if the body contour approaches a parabola. The final (tenth) 
column lists elements curvatures. 

A new output is entitled Automatic Rake Calculation. For each rake 
the two input points that bound the rake and the input identification are 
output together with three calculated quantities. The first is the surface 
area of the cone frustum defined by the tv/o input rake points. The third is 
the total flux of fluid that crosses this area per unit time. The second 
quantity, average velocity, is the ratio of the third and the first quantities. 
There is a rake output for all axisyrometric flows, both uniform and 
non-uniform, but there is no rake output for cross-flow, because its cir- 
cumferential variation guarantees zero flux. 

With regard to surface vorticity solutions, there are only two possi- 
bilities. Either there are no vorticity solutions or there is one solution 
for each body, which corresponds to a unit vorticity strength on that body 
and zero vorticity strength on all other bodies. The order of these solutions 
is the same as the order in which the bodies are input. Thus, in particular, 
an inlet with centerbody has tw*o vorticity solutions - one with vorticity on 
the inlet and one with vorticity on the centerbody. The second of these is 
not meaningful and should be discarded. (It may, of course, occur first on 
the output). 
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To help the user verify that the program is performing correctly 
with hi:-* equipment, a test case is included. The geometry is a sphere 
defined by 13 points. Axisymssetric, vortex, and crossflow solutions 
are computed using the curved element, parabolic source density, and inter- 
nally-computed curvature options. The input data and the computer output 
for this case follows. 
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PROGRAM 600F — PARABOLIC AX ISYMMFTR IC AND CROSSFLO^ 
***** CASE CONTROL DATA ***** 


12 ELEMENT SPHERE. CASE NO. 

BOOISS * 1 

NNU * 0 

CHORD * nOOQOOGO 

MACH NO.* 0*0 

TCNST * 0.0 

P$F NO. - 


SURFACE OF REVOLUTION 
CROSSFLOW 

MATRIX SOLUTION BY TR I ANGULAR 1 l AT ION (SOLVIT I 
STRIP VORTEX 

INPUT T APE* NJ* FOR COORDINATES AND NON-UNIFORM FLOW ONLY * 5 





*v 


DOUGLAS AIRCRAFT COMPANY 
LONG BEACH DIVISION 


u element sphere. 


NN » 

u 

MX * 

0.0 

MY * 

0.0 

THETA * 

0.0 

ADDX * 

0.0 

ADDY « 

0.0 

XE * 

0.0 

YE * 

0.0 




CURVED ELEMENTS! 10ECMF»QI P ! ECEWI SE-PARABOUC SOURCE DENSITIE$USIG'F*OI INTERNALLY-COMPUTED ELEMENT CURVATURES! ICURVN«C) 
NEW VELOCITY FCRMuLAE ARE USED. 

ON-SODY COORDINATES 



«**UNTRAN$FCRMED»** 

*** transformed *** 








X 

Y 

X 

Y 

X C.p. 

Y C.P. 

delta s 

SUMDS 

D ALPHA 

CURVATURE 

i 

-1.00000 

C*00000 

-1.00000 

0.00000 

-0. 99 141 

0. 13052 

0.26179 

0.26179 

0.0 

-i.OOOOX 

2 

-0,96593 

0,25882 

-0,96593 

0.25802 

-0.92385 

C. 3826 7 

0.26179 

0. 52359 

-15.00024 

-1.00000 

3 

-0.86602 

0. 50000 

-0.86602 

0.50000 

-0.79332 

0. 60874 

0.26179 

0.78538 

-14.99982 

-uooooo 

A 

-0.70711 

0.70711 

-0.70711 

0, 70711 

-0.60874 

0. 79332 

0.26175 

1.04717 

-14.99993 

-1.00000 

5 

-0.50000 

0. 86602 

-0.50000 

0.86602 

-0.38267 

0. 92385 

0.26179 

1. 30896 

-14.99989 

-1.00000 

6 

-0.25882 

0.96593 

-0.25882 

0.96593 

-0,13052 

C. 99141 

0,26179 

1.57076 

-15,00015 

-1,00000 

7 

0.0 

i. 00000 

0,0 

l.COODO 

0.12052 

C. 99141 

0.26179 

1.83255 

-14.99993 

-1.00000 

6 

0.25882 

0.96593 

0.25882 

0.96593 

0.38267 

0.92385 

0.26179 

2.09434 

-15.00003 

-1.00000 

V 

0.50000 

0, 86603 

0,50000 

0.86603 

0.60874 

0. 79332 

0.26179 

2.35614 

-14.99991 

-1,00000 

10 

0. 70711 

0. 707 U 

0,70711 

0.70711 

0.79332 

0.60874 

0.26179 

2.61793 

-15.00006 

-1.00000 

1L 

0. U6&C2 

0. 50000 

0.86602 

0.50000 

0.92385 

0.38267 

0.26179 

2.87972 

-14.99994 

-1.00000 

12 

0.96593 

0.25882 

0.96593 

0.25882 

0.99141 

0, 13052 

0.26179 

3.14152 

-15.00008 

-1.00000 

13 

1.00000 

0.0 

1.000 00 

0.0 
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12 ELEMENT SPHERE. CASE NO. 

CN-BODY UNIFORM AXISY1M6TRIC FLOW 
TRANSFORMED COORDINATES 



X 

Y 

T1 

CP 

SIN A 

COS A 

SIGMA 

N 

PHI 












£0 

1 

-1.0000000 

-0.991*083 

0. COOOQXO 
0. 1 305223 

0.1956868 

0.9617067 

0.99145 

0.13053 

-0.1183583 

-0.0000031 

0.0 

•£ & 
2 2 

> c 

2 

-0.V65926Q 

-0.9238*51 

0.2588200 

0.3826699 

0.5737293 

C.6708405 

0.92368 

0.38269 

-0.1102906 

-0.0000015 

0.0 

i 

-0.8660250 

-0.7933235 

0. 5CGOOOQ 
0.6037390 

* 0. 9128096 

0. 1667787 

0.79335 

0.60876 

-0.0947085 

-0.0000004 

0.0 

o 
t; rq 

4 

-0, 7071060 
-0.608 f 385 

0. 7071069 
U. 79 Jj^39 

• 1.1896363 

-C. 4152298 

C. 60876 

0.79335 

-0.0726728 

-0.0000004 

o 

. 

o 

* -..3 

r V *-2 

9C 

V—* 

* O 

5 

-0. 5U0G0C0 
-0. 3826694 

0.6660250 

0.9233652 

-1.3853798 

-0.9192762 

0.38269 

0.92386 

-0.0456645 

-U. 0000004 

0*0 

6 

-O.2583H0 
-0. 1305219 

0.9659260 

0.9916332 

. 1.4867077 

-1.2102995 

C. 13C53 

0.99145 

-0.0155819 

-0.0000005 

0.0 

1 

0.0 

0. ' 3052 16 

t.UOOOJJO 
0.591 6082 

1.4667067 

-1.2102966 

-C. 1 3053 

0.99145 

0.0155820 

0.0 

0. G 

»-d ^ 

o 

d 

0.2588190 
0. 3 626638 

0.96^0260 

0.9233656 

1.3853798 

-0.9192762 

-0.38268 

0.92388 

0.0436843 

-0.0000002 

0.0 

O H 
50 K 

9 

0.4*999990 
0. 6067366 

0 »►* 660260 
0. 79 3 324 D 

1.1896353 

-0.4152317 

-C. 60876 

0.79335 

0.0 726725 

0.0000002 

o 

. 

o 

Si 

10 

O.7OM069 

0.7933239 

0. 7071369 
0.6087395 

0.9128108 

0.1607765 

-0.79335 

0.60876 

0.0947084 

-0.0000001 

0.0 


u 

0. 6660250 
0.92 >8952 

0.50000 LQ 
0. 3326699 

0.5737217 

C. 6708434 

-0.92388 

0.38268 

0.1102907 

0.0000015 

0.0 


u 

0*9659260 

0.26681 9 j 






0.0000016 



0.9919082 

0. 1 33521 4 

0.1956850 

0.961 7074 

-0.9914 5 

0.13053 

0.1183582 

0. 0 


n 

i. 0000000 

0.0 










ADDED MASS 

* 0.0 

VOLUME 

« 4.1762981 

SUM (T ) t CELT A S) * 

2.9989319 
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12 ELEMENT SPHERE. CASS NO. 


ON-BCOY STRIP /ORTEX FLOW 

TRANSFORM COORDINATES 



X 

Y 

T 1 

CP 

SIN A 

COS A 

SIGMA 

N 

PHI 

l 

*1.0000000 

-0.991*083 

0*0000010 
0, 1305223 

1.335*311 

-0.7833757 

0.991*5 

0.13053 

2.2962*8* 

0.0000916 

0.0 

2 

-0,9659260 

-0,9238*51 

0 • 2 5b 3200 
0, 3326699 

-0.25*2801 

0.9353*17 

0.52338 

0,38269 

1.2829266 

0.00000*8 

0.0 

3 

-0. a660250 
-0,7933235 

0. 5 CCOGQO 
0.6C87J9Q 

-0.0182657 

0.999666* 

0.79335 

0.60876 

0.833709* 

0.0000057 

o 

• 

© 

* 

-0. 7071060 
-0.009738; 

0. 7071069 
0. 7933239 

-0.005*169 

0.9999707 

0.60876 

0.79335 

0,5359750 

0. 0000038 

o 

* 

o 

5 

-0. 5000000 
-0. 382669* 

0, 8660250 
0.9238*62 

-0.0023832 

0.99999*3 

0.38269 

0.92388 

0.3025500 

0.00000*8 

0.0 

6 

-0.2569190 
-0. 130521* 

0, 9659260 
0.591*082 

-0.0013103 

C, 9999983 

0. 13053 

0.991*5 

0.0980312 

0. 0000039 

o 

* 

o 

7 

O.U 

0. 13052 i* 

1,0000010 

0.991*082 

-0.0013151 

C. 9999983 

-0.13053 

0.991*5 

-0.0980305 

0.0000006 

0.0 

8 

0.2583190 
0. 3926o83 

0.9659260 

0,9238*56 

-0.0023851 

0.99999*3 

-0.38263 

0.92388 

-0.3025*91 

-0.0000010 

o 

. 

o 

9 

0.* 9999 90 
O.oUO 73o6 

0.8660260 

0.79332*5 

-0.005*121 

0.9999707 

-0.60876 

0.79335 

-0.5359737 

-0.0000010 

© 

* 

o 

10 

U. 7071 0e»9 
0. /S33239 

0. 7071069 
0,6037395 

-0.018271* 

0.9996662 

-C. 79335 

• 0.60876 

-0.8337085 

-0.0000029 

o 

• 

o 

U 

0.3660260 

0.9238*52 

0,5000010 
0. 1626699 

-0.25*2858 

0.9353387 

-0.92388 

0.38268 

-1.282929* 

-0.0000010 

0.0 

12 

0.9659260 

0.991*082 

0. 2 5d6 190 
0. 130521* 

1.33539*9 

-0.783279* 

-0.991*5 

0.13053 

-2.2962656 

-0.0000763 

o 

• 

o 

13 

i.UOUOOJO 

0.0 









ADDED MASS 

* 0.0 

VOLUME 

■ *.1762981 

SUM (T )( DELTA SI ■ 

0.5501668 
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DOUGLAS AIRCRAFT COMPANY 
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U ELEMENT SPHERE. 


CN-8UDY UNIFORM CRCSS FLOW 
TRANSFORMED COORDINATES 



X 

Y 

T2 

T3 

1 

-1.0000000 

0.0000010 




-0.991*083 

0. i 305223 

1 *4970138 

1. *993650 

2 


U.2W2DQ 




-0.9238*51 

0. 3626699 

1.3826275 

1. *99126* 

3 

-0. 86o02 50 

0.5000000 




-0,7933235 

0.6C87390 

L. 186730* 

1. *989996 

* 

-0. 7071060 

0. 7071069 




-0*0087365 

•0. 7933239 

0.9115025 

1. *991226 

5 

-0.0000000 

0.8660250 




-0. 3U2ou9* 

0.9230*52 

0.57*1063 

1.9993992 

& 

-J.25&8190 

U. 9059260 




-0. L3u52i* 

0.991*002 

0.1961133 

1. *996233 

7 

0.0 

t .OOOOOUO 




0. i 30521* 

0. 991*042 

•0. 19611** 

1* * 9^62 33 

8 

0.2588190 

0,9659260 




0.3826688 

0.9233*56 

-0.57*1080 

1. *993992 

9 

0. *'9*59590 

0. 86OU260 




0.608 73d 6 

0.79332*5 

-0,9115077 

1 . *99 1217 

10 

0. 7071069 

0, 7 07 10o9 




0. 7933239 

0.6087395 

-1.1867332 

1. *989956 

11 

0. 86602 ^0 

0. 5000010 




0.9238*152 

0. 3326699 

-1.3826275 

1. *991226 

12 

0.^659260 

0.258dl90 




0.991*082 

0. 1 3J521* 

-l. *070119 

1. *993696 

13 

1.0030000 

0.0 




CASE NO. P$F « 


SIN A 

COS A 

SIGMA 

N 

p,Hl 

0.991*5 

0.13053 

0,0155821 

0.000000* 

0, 065178* 

0. 92388 

0.38269 

0.0*5605* 

0.0000003 

0. 1910008 

0. 79335 

0.60876 

0.072*52* 

0.0000006 

0.3037606 

C. 60876 

0.79335 

0. 09***96 

Q.QG0QQ03 

0.3959663 

C. 38269 

0.92368 

0, 1101502 

-0.0000001 

0. *613678 

C. 13053 

0.991*5 

0.1183575 

0.0 

0. *953313 

0,13053 

0.991*5 

0.1184576 

0.0000001 

0, *933312 

•0. 38268 

0.92388 

0.1101502 

O.OCOOOOl 

0.4613683 

C. 60076 

0.79335 

0. 09***93 

0,0000002 

0.3959657 

■0. 79335 

0.60876 

0.072*517 

0.0000003 

0.3037609 

C. 92388 

0,33268 

0.0*560*7 

0.0 

0. 190999* 

0. 99 L*5 

0. 13053 

0.0155822 

0.0000002 

0.0651785 


ADDEC ’ MASS - 


2.08560*7 


VOLUME • 


*•1762981 


reproducibility op the 

ORIGINAL PAGE IS POOR 
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APPENDIX: 


USAGE OF THE IPTANV FLAG 


Usage of the IPTANV flag permits tangential velocity to be specified 
over part of the surface of the body or bodies about which flow is being com- 
puted. On the remainder of the surface the normal velocity is made zero in 
the usual way. At locations where tangential velocity is specified the normal 
velocity will not in general be zero. This feature can be used only for axi- 
symmetric flows with a uniform free stream. 

If the flow about several bodies is to be calculated, these bodies are 
considered in sequence as determined by the order in which they are input. 

This feature uses a flag to legislate that normal velocity conditions will 
be applied on the first several bodies of this sequence and tangential ve- 
locity conditions will be applied on the remaining bodies. All bodies on 
which normal velocities are specified must be input before any of the bodies 
on which tangential velocities are specified. For example, suppose it is de- 
sired to specify tangential velocity over the region near the middle of a 
body, while the usual normal velocity conditions is applied over the forward 
and aft regions. This body mus** be input as three bodies, the first two 
cf which are the forward and aft regions and the third of which is the middle 
region. 

The desired tangential velocities are input as a set of numbers. The 
first number in the set is the tangential Velocity on the first element of 
the first body on which tangential velocity is specified. The subsequent 
numbers of the set give tangential velocitites on the remaining elements in 
the order that they were input. The order of the input of the tangential 
velocities is used to associate them with elements, and it is the responsi- 
bility of the user to insure that this order is correct and that the number 
of input tangential velocities equals the number of elements ca the bodies 
on which tangential velocities are specified. The tangential velocities 
should be input as total tangential velocities, including the contribution 
of the uniform free stream. 

If the desired tangential velocity is a constant over the body or bodies 
on which tangential velocity is specified, only this single constant need be 
input. The constant should represent a total tangential velocity. It is not 
possible to input a separate constant for each body on which tangential 
velocity is specified. 

It is possible to specify tangential velocity over the entire surface 
of the body or bodies that are input, so that no normal velocity condition is 
applied. 

The final output is similar to that of the usual case. On the bodies 
where tangential velocities were specified the printed values of T should be 
the same as those input. The normal velocities will not in general be zero, 
but their magnitude is a measure of how near the input shape is to an impervious 
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body having the same tangential velocity distribution. The printed values of 
N are outward normal velocities. A plus sign denotes flow out of the body* 
and a minus sign denotes flow into the body. 

The option of specifying tangential velocity is excercised by punching 
some non-zero integer into C.C. 14 of the control flags card (IPTANV). 

This value indicates the number of bodies on which tangential velocity is 
being specified and therefore must be less than or equal to the total number 
of bodies being input. In particular, if a three body case is being run 
and tangential velocity is to be specified on two of the bodies, then those 
two bodies must be input last in the input sequence, and IPTANV must have a 
value of 2. The use of IPTANV affects the matrix solution portion of the 
program as follows : 

Whereas the usual equations to be solved for a for uniform axi— symmetric 
flow are 

N 

A.. a. = sin a (for i = 1, 2 N) 

3 i 

j +1 

Where N is the total number of elements on all the bodies, use of 
IPTANV directs the program to solve for a in the following ways: 

I. If IPTANV is equal to the total number of bodies input, the 
equations to be solved for a become 


B ij a j 23 TG i = cos a ± (for i - 1, 2, N) 

where B^j is the usual axi-symmetric Bfj matrix and TG^ is an 
array of N numbers consisting of either (i) the same constant 
number (this number to be specified by the user in the TCNST 
field of input card no. 3), or (ii) N values input from cards 
(see card set 12) . 

II. If IPTANV is less than the total number of bodies (say, IPTANV 

= k), then tangential velocities are to be specified on the last 
k bodies input. That is, if M is the total number of elements on 
the last k bodies, then N-M is the number of elements on the first 
NB-k bodies (where NB is the total number of bodies Input) , and 
the equations to be solved for a become 
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N 



j + 1 



/ 



N 

A i j a j = sin (for i = 1, 2, •••> N-M) 

j- 1 

^ ^ B i j aj = TG i - cosa. (for i - N-K+l, . . . , N) 

If all values of the TG array are input by cards, the user must 
take care k to input precisely H values. 

i 

Note that his special matrix formation finally affects only the 
values of a for axi- symmetric cases (crossflow matrix formation 
and solution remains unchanged) ; all remaining calculations are 
performed as usual. 

Restrictions exist in the use of IPTANV. In particular, IPTANV 
cannot be used in conjuction with any non-uniform flow solution, 
whether axi-symmetric or crossflow. Neither can IPTANV be used 
if the total number of body elements exceeds 275. Nor can the 
strip vortex option (IVORT) be used. 
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